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We present the theoretical basis for and experimental verification of arbitrary single-qubit 
state generation, using the polarization of photons generated via spontaneous paramet- 
ric downconversion. Our precision measurement and state reconstruction system has the 
capability to distinguish over 3 million states, all of which can be reproducibly gener- 
ated using our state creation apparatus. In order to complete the triumvirate of single 
qubit control, there must be a way to not only manipulate single qubits after creation 
and before measurement, but a way to characterize the manipulations themselves. We 
present a general representation of arbitrary processes, and experimental techniques for 
generating a variety of single qubit manipulations, including unitary, decohering, and 
(partially) polarizing operations. 

Keywords: photon polarization, qubit, quantum process, tomography, parametric down- 
conversion 

Communicated by: J Shapiro & H-K Lo 

1 Introduction 

Quantum information processing (QIP) promises great potential power over its classical coun- 
terpart in the areas of computing, communication, and metrology jj. Nearly all QIP protocols 
require specific initial states and the ability to manipulate the qubits with exquisite precision. 
Here we demonstrate the creation of arbitrary single qubit states encoded in the polarization 
of single photons. Single-photon Fock states are conditionally realized by detecting one photon 
of a pair produced in the process of spontaneous parametric downconversion [2] . We manip- 
ulate the photons' polarization state using a series of birefringent waveplates (which enable 
any unitary transformation) and a thick birefringent decoherer (which allows the production 
of mixed states). Using a method of state tomography, we can experimentally determine the 
most likely density matrix which describes the resulting quantum state. Finally, we have 
implemented several versions of process tomography, by which we can accurately characterize 
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any quantum process acting on the polarization qubit, including all unitary transformations, 
(partial) measurements, and decoherence. 

2 State Creation 

We represent our qubits by single-photon polarization states in the horizontal- vertical basis, 
with the state vectors 

|0) = \H) = ( J \ and |1> = \V) ee ( J ) , (1) 
or in density matrix notation, 

PH = ( I I ) and Pv = ( I I ) " (2) 

The density matrix of an arbitrary single qubit can be represented by three independent real 
parameters (A, B, and 6): 

P= Be -is X _ A , (3) 



where 0< A <1 from normalization, and \B\ < y A(l — A) from positive semi-definiteness. 
Another equivalent representation is given by 

P=\(I + r-a), (4) 
where I is the identity and we define the polarization analogs of the Pauli spin matrices as 

ai = ( i I ) ' ff2 - ( i o J ' and a3 = ( -1 ) ■ (5) 

The components of r give the degree of polarization for the photon in the Horizontal- 
Vertical (H-V), Diagonal- Antidiagonal (D-A), and Right-Left Circular (R-L) bases? As such, 
they are often identified as rjj = T\, rjj = ra, = r 3 ? It is useful to view these components 
as coordinates in a 3-D space of polarizations; the constraint |r| < 1 implies that all states 
must lie inside or on a sphere of unit radius, known as the Poincare sphere. Points on the 
surface of the sphere (|r| = 1) represent pure polarization states (linear polarization states on 
the equator, right and left circular polarization on the north and south pole, respectively), 
while points inside the surface (|r| < 1) represent partially mixed states. The center of the 
sphere (|r| = 0) corresponds to a completely mixed state, i.e., an unpolarized photon. 

Because an arbitrary state has three independent parameters, the generation of such a 
state requires three adjustable elements. Considering the Poincare sphere, we make an ansatz 
that a half-waveplate (HWP), followed by a thick birefringent decoherer, a half-waveplate, 
and a quarter-waveplate (QWP) are sufficient to generate all one-qubit polarization states 
from a pure linear polarization fiducial state (\H) in our case). We now derive formulae that 



"These components ri are related to the Stokes parameters (Si) of classical optics IH1 bv n = j-j, 

6 The conversion between the representations in 131 and l%l is given by r^j = 2A — 1, rrj = 2Bcos(<5), and 

r R = 2Bsin(5). 
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give waveplate settings for an arbitrary state, thus proving our ansatz. The operators that 
represent the action of half- and quarter-waveplates, respectively, are the Jones matrices 0] 

_ / -cos26» -sin26> \ 

O hwp (6)=^ _ s . n2e cQs2e j (6) 

and 

„ lCi \_( 1 - (1 + i) cos 2 9 — (1 + i) sin cos \ . . 

O qwp {9) = ^ _ (1 + z)sin0cos0 l-(l + i ) s in 2 )■■ (7) 

where in each case the parameter 9 is the angle the optic axis makes with horizontal. 

To create arbitrary states as in 0, we start with photons in the state P h, and direct them 
through a half- waveplate at 9i, giving 

-r» ia\ n to \\ ( cos 2 29 1 cos 29 1 sin 29 1 \ . , 
Pl = Owp^O^O^p^Jt = cos2 ^ sin20i sin 2 20i ) = |2(9i)(2(9 1 |, (8) 

i.e., the pure linear polarization state \29i). As shown in the first box in Fig.Q this operation is 
described on the Poincare sphere by rotating the state \H) by 180° about an axis - representing 
the optic axis of the waveplate - that lies 29 1 away on the equator; the factor of 2 arises because 
\V) lies on the opposite side of the Poincare sphere, i.e., 180° away from \H)9 

The next step is to introduce decoherence by separating the horizontal and vertical po- 
larization components by an optical path length difference much longer than the coherence 
length of the photons (see Appendix A) 0. If the coherence length is longer than a few 
millimeters, one can use a polarization-dependent delay line, as shown in the second box of 
Fig. n For our downconversion source, interference filters at the detector typically define the 
spectral bandwidth; for filters with 10-nm full width at half maximum, the coherence length 
is only 49-/^m. In this case it suffices to use a thick birefringent element [e.g., 1 cm of quartz] 
to completely decohere the polarizations within the eigenbasis of the decohering element UJ. 
We can control the amount of decoherence by tuning the polarization of the input state. 
For example, if the polarization before entering the decoherer is \H) (or |V)) then the state 
purity is preserved; if the input state is diagonally polarized, i.e., \D), the resulting state is 
completely mixed. An arbitrary value of |r| is produced by setting the orientation of the first 
half- waveplate to 9\ = | arccos |r|. After the rotated light is directed through the birefringent 
decohering element, the reduced density matrix describing only the polarization is of the form 
(see Appendix A): 

, / cos 2 29 1 \ 
Pi=[ sin 2 29, )■ (9) 

Next, using waveplates, we unitarily transform JUJl into the desired final state. Note that 
iJSJ can be rewritten as an incoherent sum of a horizontally polarized pure state and an 
unpolarized, completely mixed state: 

p\ =cosA9 l \H){H\+2sm 2 29 lPmlxedl (10) 

where p m ixed = 5I is the completely mixed state. Because quantum mechanics is linear, we 
may operate individually on each part of this sum to realize the final state. As the unpolarized 

C A very complete discussion of the use of the Poincare sphere to describe the action of various crystal optics 
may be found in 
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Fig. 1. The experimental setup to realize an arbitrary single (polarization) qubit, along with the 
representation of the state in the Poincare sphere at each step of state preparation. In the first 
stage, on the left, horizontally polarized photons are sent through a half-waveplate (HWP) with 
an optic axis at 9\ from horizontal, giving the linear polarization state |20i). This is then sent 
through a decoherer that has an optical path length difference greater than the coherence length of 
the photon. One method to achieve such a decoherer is with a polarizing beamsplitter that sends 
vertical polarization in a delay loop while horizontal is transmitted, thereby suppressing the phase 
coherence in the horizontal- vertical basis (see Appendix A), and eliminating any amplitude in the 
off-diagonal elements of the state. Decoherence in the Poincare sphere appears as a projection 
of the state onto the \H) — \V) "spindle" (second box). By appropriately adjusting the first 
HWP, states with arbitrary mixedness can be produced, ranging from a pure state (|r| = 1) to 
a completely mixed state, i.e., unpolarized photons. The last two waveplates, a HWP and a 
quarter- waveplate (QWP), set the final direction of the (possibly mixed) state on the sphere. 



part is is unchanged by any unitary transformation, the final form of the state is determined 
by operating on the \H){H\ term with a half-waveplate (at 62) and a quarter-waveplate (at 
6* 3 ). Algebraically determining the desired values of 62 and 83 to obtain a particular final state 
is non-trivial. However, by noting the geometric action of these waveplates on the Poincare 
sphere (see the third and fourth boxes of Figure 0), one can arrive at the following final 
solutions for the waveplate angles* 



5[V(2A-1) 2 +4B 2 ] 



(11) 



d In terms of the Poincare sphere vector r, we have 



03 



— arccos (|r|) 



arctan ( —2- ) + arctan ( - ' 
r H 



— arctan ( — — ) . 
2 r H 



e In any real system, the various waveplates may have slight errors in their retardations. Nevertheless, we have 
shown that any single qubit state may be produced, as long as the QWP retardance error is smaller than the 
HWP retardance error. In this case, however, we must numerically search for the optimal waveplate settings. 
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Fig. 2. Single-qubit downconversion experiment. Detection of a photon in the upper arm condi- 
tionally prepares a single-photon state in the lower arm. The polarizing beamsplitter (PBS) after 
the nonlinear BBO crystal prepares an extremely pure initial state \H), which is then transformed 
by the waveplates and decoherer (~1 cm quartz birefringent element). The tomography system 
allows an accurate measurement of the resulting density matrix p (see Sect. 3). 
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Armed with these equations we set up our experiment as in Figure |3 The process of 
spontaneous parametric downconversion conditionally prepares single-photon input states ■ 
A nonlinear crystal (BBO) is pumped with a vertically polarized 351-nm beam from an 
Argon- ion laser (average power 86 mW). The BBO is cut to produce non-collinear frequency- 
degenerate horizontally polarized photon pairs at 702 nm via type-I phase matching. The 
non-collinear photon pairs are collected in two modes, separated by 6° outside of the crystal. 
The first mode impinges onto a detector assembly consisting of a 10-nm FWHM interference 
filter (IF) centered at 702.2 nm, a lens, and an avalanche photodiode (APD) operated in 
Geiger mode (Perkin-Elmer #SPCM-AQR-14). Single-photon Fock states are prepared in 
the second arm conditional on a "trigger" count in the first detector 0. The second mode 
passes through a polarizing beamsplitter (PBS) (to ensure a good fiducial state \H)), the state 
preparation waveplates-^and decoherer (1 cm quartz slab) discussed above, a QWP-HWP-PBS 
combination for state analysis (see Sect. 3), a 2.2-mm iris, and a detector assembly identical 
to that in the first mode. 

We created a variety of single qubit states using this system. Some sample reconstructed 
density matrices are shown in Figure [3J One measure of our ability to accurately prepare 
specific states is the fidelity 9 between the measured and the target states. We typically 
observed fidelities better than 0.997; this number is presently limited by counting fluctuations 
due to Poisson statistics, and also to small drifts (less than 0.5%) in either the pump laser 
intensity or the detector efficiencies. For the data in Fig. |3| a total ~150,000 counts were 
accumulated for each state. 



^All of the waveplates are in motor-controlled stages, and can be set with an accuracy of less than 0.1°. 
ff The fidelity is a measure of state overlap 0. For two pure states, \ipi) and \ip2), the fidelity is simply 
FQifii), \ip2)) = [v^lhfei)! 3 ! while for two general density matrices, p\ and p%, the fidelity is -F(pi,P2) = 
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Fig. 3. Single qubit data. Shown are the target density matrices, plots of the target and experi- 
mentally measured density matrices, and the fidelity (F) between them. The imaginary elements 
of the density matrices are not shown because they have zero amplitude for the target states, and 
are always less than 2% for the experimental states. 



3 State Characterization 

One cannot determine the unknown polarization state of a photon with a single measurement. 
Instead, a large ensemble of photons prepared in an identical manner must be projected into 
different polarization basis states. Four measurements are needed - three to determine the 
relative values of the three independent parameters that characterize an unknown state, and 
a fourth to determine normalization. The polarization analysis is carried out setting the 
quarter- and half-waveplates shown in the tomography system box of Fig. Although there 
are many choices as to the particular measurements that can be made (in principle measuring 
along any three non-collinear bases suffices), we choose the analysis states: (V>o| = (V\, 
(-01 1 = (H\, (^2 1 = (D\, and (-03 1 = (R\ corresponding to N , N\, N 2 , and N 3 photon counts 
in some fixed measurement time interval (typically 100 s). As explained earlier, to obtain 
single photon Fock states, we count in coincidence; thus, JV, are coincidence counts between 
the two detectors shown in Fig. [21 Since we need to estimate probabilities, we must measure 
a complete basis to normalize the photon counts. As \H) and | V) form a basis, Af = No + N% 
gives a convenient normalization factor for the r^. r# = 2N\/M — 1, td = 2N2/N — 1, and 
Tn = 2N3/N — 1. The density matrix of the state may then be reconstructed as in J2J. 

Unfortunately, as this state reconstruction is based on photon counting, statistical fluctu- 
ations or drift often yield an unphysical result; therefore, we employ a maximum likelihood 
technique to estimate a physical density matrix that would most likely produce the measured 
data. James et al. describe the technique for determining the joint state of two qubits [S]. Here 
we distill their argument down to the one qubit case (if states other than single-photon states 
are used, more sophisticated methods must be used to determine the quantum polarization 
state 0). 

A physical density matrix representation has three cardinal properties: normalization, 
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positive semidefiniteness, and hermiticity. Therefore, we first guess a density matrix that by 
definition has the aforementioned properties. A matrix that can be written as T^T is positive 
semidefinite and hermitian 8 . To normalize such a matrix, we divide by its trace so that 
T^T/Tr(TtT) has the three properties for a legitimate physical density matrix. We choose 
the following invertible form for T(t): 

T(t) S T(t 1 ,* a ,* 8 ,t 4 ) S ( ts * 1 . t4 I). (14) 

Using 114(1 . the density matrix formula is 

P P {hMMM) = T+(t)T(t)/Tr(Tt(t)T(t)). (15) 

Next we will introduce a likelihood function that quantifies how similar p p (tx, t2, £3, £4) is to 
our experimental data: 

C(N N N N-t ft f) _^ W(j>M^ 2 ,t 3 ,t A Mu)-N v ] 2 

^ 2A/(Vv|/M*i>*2,£3,£4)|Vw 

where the quantity Af{if>v\pp(ti, t 2 , £3, t&)\ij) v ) represents the expected number of counts for a 
projection of our trial density matrix p p onto the analysis state |Vv)- Note that the coincidence 
counts are subtracted from J\f(ip v \p p (ti, t 2l ti,ti)\ip v ) so that the likelihood function must be 
minimized to obtain the set of t^° pt ^ , and therefore the state that best retrodicts the actual 
measurements^ 

The question of how many different states can be reliably produced is limited by how 
many may be reliably experimentally distinguished 1 . To determine the number of accessible 
distinguishable states in our system, several states in different parts of the Poincare sphere 
(with approximate |r| values of 0, 0.25, 0.5, 0.75, and 1) were created and measured 10 
times each. We measured for 10 minutes per tomography yielding 300,000 counts per state 
reconstruction. For each set of measurements the average state vector r was determined. 
Next, the standard deviation of the ten trials was determined for three directions: f, and 
two directions transverse to f. Assuming a Gaussian distribution along these directions, 
we calculate that an ellipsoid with semi- axes equal to 1.69 times these standard deviations 
will account for over 95% of the events. These uncertainty ellipses are mainly due to count 
fluctuations from laser and detector efficiency drift, as well as intrinsic Poissonian counting 
statistics. One unexpected result of our measurements was that the thickness of the ellipsoid, 
i.e., the length of the minor axis along the direction of r, depended on the radial coordinate 
|r|, and varied from a minimum value of 0.0021 for |r| = 1, to a maximum value of 0.0062 for 

'"The minimization is performed using the MATLAB R12 function fminunc which requires an initial estimate for 
the ti. For the initial estimate, starting guesses for the pi are calculated from the measured counts No, Ni, N2, 
and N3; through ||}, this gives a starting estimate for p, which, from 1141 and 1151 . yields starting values for 

the ti : tl = 1 - v /l- 5JI i- 7 ( ( l- rH )2 + r 2 )+r 2 ) , t2 = VWElH , t3 = , and u = 

l When specifying a produced quantum state, it is an interesting question whether one should give an error 
ball (e.g., a patch on the surface of the sphere for a pure state), or simply average over this ball to yield 
a slightly mixed state. The former approach accounts for the fact that the error ball could shrink if more 
data were taken. The latter method acknowledges that the density matrix already encodes the totality of our 
knowledge, based on our measurements on identically prepared members of an ensemble. 
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Fig. 4. Filling the Poincare sphere with experimentally determined uncertainty "patches" . These 
arise primarily due to counting fluctuations, either from the Poisson statistical uncertainties in 
photon counting, or from slight drifts (< 0.5%) in the laser power or detector efficiencies. The 
measured uncertainty patches are approximately ellipsoids (pancakes) as shown above; for clarity 
their dimensions are shown scaled up by a factor of 5. Note that the thickness of the pancakes 
depends somewhat on the mixture of the state, indicating that some regions of the Poincare sphere 
are more sensitive to counting statistics. The size of uncertainty patches implies that we are able 
to distinguish more than three million unique single-qubit states, assuming ~10 minutes data 
collection time per state, i.e., ~300,000 detection events). 

r| = 0.25 (see Figure0J. Numerical simulations support this observation. Taking into account 
the varying size of the uncertainty patches, and assuming an approximate close packing of 
the entire Poincare sphere volume, we estimate that we can reliably distinguish over three 
million states. 

4 State Manipulation 

In order to gain total control over a single qubit, it is necessary to have an understanding of 
what operations are physically possible. Any physical process acts as a map from all possible 
input states to a transformed set of output states. For example, the identity process maps 
every possible input state to itself. In other words, it has no effect on any input states. Other 
familiar processes include the unitary transformation, which when viewed in the geometric 
picture of the Poincare sphere, is simply a rotation around a fixed axis (Fig. [5jt) ; projections 
which cither partially or completely project a state into a single basis (a single point on the 
Poincare sphere, see Fig.^Ja); or (partial) decoherence, which partially or completely collapses 
all states to the "spindle" formed by the eigenstates of the decohering interaction (Fig. [5^). 

Note that while some processes, such as the identity, projection, or the unitaries, can 
be represented by a single 2 by 2 matrix (which operates on a two element column vector 
representing the qubit, as in Jones calculus), some more complex processes such as decoherence 
or incoherent partial polarization require a different method of characterization. Consider a 
simple operation (e.g., unitary or projecting) E acting on a general state pi„. The output 
state is given by 

Pout = Ep m E J< . (17) 
A general operation £ cannot be represented using a single E matrix, but instead can be 
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Fig. 5. Shown are several example single-qubit quantum processes. Each process is represented 
by both a geometric picture and by operator matrices. The geometric picture shows initial states 
(white balls) on the Poincare sphere mapped to final states (solid balls). The effect of a process 
£ on an input state can be defined as £(pi„) = Ejpj n E\. The operator matrices shown 
correspond to the Ej matrices, (a) An example of a unitary transform. This transform, the 
Hadamard gate, corresponds to a 180° rotation about the 22.5° axis, (b) A total projection 
operator. This corresponds to a perfect horizontal polarizer, (c) A coherent partial polarizer. 
This partial projection corresponds to a device which transmits all horizontal light and half of 
vertical light. These two components maintain their phase relationship. Note that pure input 
states remain pure while all states travel towards H on the sphere, (d) An incoherent partial 
polarizer. This corresponds to an operation which 50% of the time projects into the H basis 
and 50% of the time does nothing. The photons which are subjected to the first operation are 
incoherently added to those of the second operation; thus this process transforms pure input 
states into mixed states, (e) A decoherer. This operation decoherers in the H-V basis, removing 
the phase relationship between these two components. All states travel in a straight line to the 
spindle connecting the H and V states. 
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described using the following operator-sum decomposition pQ 

Pout = £{pin) = ^EjPinEj- (18) 

j 

An arbitrary number of E operators can be used to represent a process, and as we saw above, 
decohering and partially polarizing processes require at least two matrices. The necessity of 
multiple E operators hints at the variety of processes which must be possible, and the subtlety 
involved in representing an arbitrary process. Consider, for example, the partial polarizers 
shown in figures |SJ: and [SJi. Figure [5]; shows a coherent partial polarizer, which acts as a 
partial projector. Only a single E matrix is necessary to represent it. Experimentally, this 
projector is implemented using a series of tilted glass plates for which Ty /Th is 0.5-? All of 
the vertical light which passes through the glass plates maintains a definite phase relationship 
with all of the horizontal light which passes through. For this reason we call this a coherent 
partial polarizer. 

Now consider a process which totally projects into the horizontal basis, but only acts 
on 50% of the measured qubits (photons). The light passing through the polarizer has no 
phase relationship with the light that does not, and an incoherent mixture results. We call 
this process an incoherent partial polarizer (Fig. 01). In fact, any process which transforms a 
pure state into a mixed state must be represented by having one E operator acting on a certain 
percentage of qubits, and a separate operator or operators acting on the remaining qubits. 
Decoherence acts in exactly this way. For example, consider the case in which the input qubits 
are randomly and incoherently subjected to either a a z rotation (in the H-V basis) or the 
identity matrix. Combining the photons from the first process with those from the second, 
we see that any coherence between the H and V bases has been destroyed, collapsing output 
states to the H-V spindle on the Poincare sphere. This action is illustrated in Fig. [SJi. 

Using equation (|18|l . we can simplify the general representation of a process. First consider 
that the E operators above can be expanded into a linearly independent basis of 2x2 operators, 
such as the identity matrix (<7o) and the Pauli matrices (<7i, (72, 03): 

3 

E^Y^CiVi- (19) 

Substituting into equation l|18|l and combining terms, 

3 

£(Pm) = ^ O-ipinCTjXij- (20) 

The Xij matrix a 4x4 positive Hermitian complex number matrix, shares many sim- 
ilarities with a state density matrix p. (For an earlier, alternate characterization of quantum 
processes, see ^J). While a density matrix shows the coherent and incoherent combinations 
of four orthogonal state vectors that make up a state, the \ matrix shows how an arbitrary 
process is represented by coherent and incoherent combinations of four 2x2 operators (e.g., 
the a matrices). In the case of the density matrix, a basis change can yield the same state 

J When the plates are tilted so that the photons are incident at Brewster's angle 56° [3J, the p-polarization is 
completely transmitted while the s-polarization suffers a reflective loss of 15% per air-glass interface. 
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written in a different orthogonal four-element basis. There is an exact analog to this orthog- 
onality condition for the \ matrix, with orthogonality defined as Tr{E\E2) = 0. In fact, if 
the x matrix is diagonalized, it becomes clear that an arbitrary process can be represented by 
the operator-sum decomposition above, with only four orthogonal 2x2 operators in the sum. 
Note that in general these four operators will not simply be the Pauli matrices o~j. 

A physical picture of this representation corresponds to four different orthogonal 2x2 
operators, each of which has a specific probability to act on an input qubit, i.e., they are 
applied incoherently. In this way an ensemble of identical input pure states can be transformed 
into a mixed output ensemble (see Fig. |SJ). 




Fig. 6. An ensemble of single qubits are operated upon by an example of an arbitrary process. In 
general, any process can be represented by four orthogonal 2x2 operators, each of which has some 
percentage chance of affecting an input qubit. These four operators have no restrictions, but they 
must in some way incoherently act upon input qubits. In this way input pure states affected by 
more than one operator can be transformed into a mixed ensemble. 



Now that we have a useful representation of arbitrary processes and a physical interpre- 
tation of that representation, it is necessary to be able to measure and reconstruct a given 
unknown process. The measurement of an unknown process is accomplished using quantum 
process tomography. The simplest way to measure a process involves preparing a complete 
basis of four single-qubit input states (e.g., H, V, D, R), subjecting each of these to the un- 
known process, and measuring the four output states. This is referred to as standard quantum 
process tomography (SQPT) ^J^^- While this is effective, it requires four input states. It 
is possible, however, with a single 2-qubit input state to exactly characterize an unknown pro- 
cess [121 1131 By using a second ancillary qubit, highly correlated with the primary qubit 
(which is subjected to the unknown process), measurements taken in coincidence on the out- 
put 2-qubit state allow reconstruction of the single-qubit process (this method is referred to 
as ancilla-assisted process tomography - AAPT). Because this technique requires correlation 
between the primary and ancillary qubits, maximally entangled states yield the most accu- 
rate AAPT results (this special case of AAPT is referred to as entanglement-assisted process 
tomography - EAPT). Surprisingly, there is a class of separable - completely unentangled - 
states which possess the necessary correlations to perform AAPT |X 51 fTO| . All three types of 



12 Precise creation, characterization, and manipulation of single optical qubits 



QPT have been experimentally realized: SQPT, first in NMR |17l IT%| and later in photon 
systems [H]; EAPT [231 ED; and AAPT using unentangled states [Ej. Shown in Fig. are 
experimental arrangements for these three techniques as demonstrated in 




Fig. 7. Experimental setups for quantum process tomography. A, B, and C above denote 
which elements are present for SQPT, EAPT, and unentangled AAPT, respectively, (a) Standard 
quantum process tomography: Polarizer (P), half waveplate (HWP) and quarter waveplate (QWP) 
allow preparation of required pure single-photon states; identical elements allow tomography of 
the post-process states, (b) Entanglement-assisted process tomography: The source produces the 
maximally entangled state (\HH) — \VV)) / A two-photon tomography of the output allows 
reconstruction of the process, (c) Ancilla-assisted process tomography: The source produces a 
separable, or completely unentangled, input state. Correlations in this state still allow AAPT. 

Experimentally, we have realized a variety of processes, allowing us considerable freedom 
to manipulate our arbitrarily generated single-qubit states. Arbitrary unitary transformations 
can be generated with a series of birefringent waveplates, specifically, two quarter waveplates 
and a half waveplate |22| . Incoherent full or partial projection is accomplished by using polar- 
izers, which are inserted for a fraction of the total measurement time commensurate with the 
strength of the desired projection. Coherent partial projection is generated using tilted glass 
plates. Perhaps the most interesting process, given current interest in quantum computation, 
is decoherence. As described in Sec. 2 and Appendix A, we introduce decoherence by pass- 
ing photons through thick pieces of birefringent quartz. By varying the thickness of quartz, 
an arbitrary strength of decoherence can be introduced. By adding unitary transformations 
(waveplates) before and after the decoherer, this decoherence can be applied in any basis. 

Experimentally measured processes corresponding to the theoretical examples given before 
are shown in Fig. [5] As before, the geometric picture offers a convenient way to visualize an 
otherwise complex operation. In fact, considering that an arbitrary process includes 15 inde- 
pendent parameters (12 to define four orthogonal 2x2 operators and 3 to assign probabilities 
to them), the ability to visualize its operation at all is somewhat surprising. 

The natural extension of this work is toward control of multiple qubits. As this work 
shows, in these exponentially larger Hilbert spaces it will be increasingly important to find 
geometric or otherwise intuitive methods to solve the much harder problems of multiple qubit 
creation, characterization, and manipulation. 
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Fig. 8. Examples of experimentally generated and characterized processes are shown. The 
dark colored mesh corresponds to the output mapping of all input pure states. The small solid 
spheres correspond to the cardinal points: H (green), V(yellow), D (purple), A(blue), R(red), and 
L(pink). The 2x2 operator matrices correspond to the Ej matrices which have a greater than 
1% contribution to the measured process. Note that in contrast to the operators in Fig. 6, 
here we have included the probability weightings directly in the matrices shown, (a) A unitary 
transformation implemented using a birefringent waveplate. (b) A horizontal polarizer. Note 
that a single solid cardinal sphere is not mapped to H. This sphere corresponds to the vertical 
input case. Some small amount of vertical polarization leaks through this measured horizontal 
polarizer. Due to the extremely low number of vertical photons which survive this process, the 
experimentally calculated output state is essentially random. This highlights one disadvantage in 
using this geometric picture: intensity information is not directly visible, (c) A coherent partial 
polarizer. Implemented using tilted glass slides, this polarizer maintains the coherence of output 
states: thus, pure states remain pure, but slide along the surface of the sphere towards H. (d) 
An incoherent partial polarizer. This was simulated by inserting a horizontal polarizer into the 
beam 50% of the time, (e) A decoherer. A thick piece of birefringent quartz causes the H and 
V polarizations to separate from each other in time. They then become partially distinguishable, 
and pure superpositions of H and V become mixed. 
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Appendix A 

Decoherence always arises from the entanglement of the quantum system being considered 
to some other quantum degree of freedom, which is then traced over. In our system we 
realize decoherence by coupling the frequency and polarization degrees of freedom and then 
measuring in a frequency- insensitive way [Sj. The state of a photon written in terms of its 
polarization and frequency spectrum is |£) = (a\H) + f3\V}) ® f cLuA(u))\ll) >, where a and 
j3 are complex normalized coefficients, and A(lo) is the complex amplitude for frequency to, 
normalized such that J duj\A(uj)\ 2 = 1. 
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To decohere in the |.ff)/|V) basis, we send the photon through a birefringent element 
whose fast axis is parallel to the horizontal polarization. Traversing a birefringent element 
of length L adds a phase that is polarization and frequency dependent, producing the state 
|£d) = J duj{e l c u a\H) + e l c ^ (3\V)) A(u)\uj >. Tracing over the frequency gives 



and the polarization state will be fully decohered. Note that although this form of decoher- 
ence due to dcphasing is reversible (i.e., by using a compensating birefringent element), it is 
fundamentially no different than any other type of decoherence, which in principle is always 
reversible if one could access the entire Hilbert space describing all parts of the experiment. 



fc The coherence length (L c ) of the photon is determined by the Fourier transform of the spectrum A(u>). For 
example, if A(u)) is a Gaussian of full width at half maximum Aw, then L c = ■ 




(A.l) 



where 4>(lj) = (nn — nv)-w. As long as (uh — ny)L is much greater than the photon's co- 
herence length^ the off-diagonal elements of expression l|A.l(l will effectively average to zero, 



